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STATIONARY DISTRIBUTION FOR DIOECIOUS BRANCHING 
PARTICLE SYSTEMS WITH RAPID STIRRING 

By Feng Yu 1 

University of Oxford 

We study dioecious (i.e., two-sex) branching particle system mod- 
els, where there are two types of particles, modeling the male and 
female populations, and where birth of new particles requires the 
presence of both male and female particles. We show that station- 
ary distributions of various dioecious branching particle models are 
nontrivial under certain conditions, for example, when there is suffi- 
ciently fast stirring. 

1. The particle models. We consider a type of particle system that can 
be used to model sexual reproduction of a certain species. This work was 
inspired in part by Dawson and Perkins [1], which studied the following 
system of stochastic partial differential equations: 

-£(t,x) = -Au(t,x) + (-fu(t,x)v(t,x)) 1/2 Wi(t,x), 

-£(t,x) = -Av(t,x) + ( 1U (t,x)v(t,x)) 1/2 W 2 (t,x), 

where A = J2id 2 /dx 2 is the Laplacian, 7 > and Wi(t,x) (i = 1,2) are 
independent space-time white noises on M + x M. One can associate u(x,t) 
and v(x, t) with the male and female populations of particles, respectively, at 
spatial location x and time t. Loosely speaking, the above SPDE system says 
that individual male or female particles move around according to Brownian 
motion, but branching is only possible when both male and female particles 
are present at the same spatial location. Notice that at spatial locations 
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where the female population is 0, the branching rate for the male population 
is also 0, therefore the male population does not die and the only effect 
on the male population at those spatial locations is the diffusive effect of 
the heat kernel. This behavior is not very realistic, since one would expect 
a natural death rate for male particles, even without the presence of any 
female particles. In this work, we study models involving a finite number of 
male and female particles with more realistic behavior. 

The model we study involves two types of particles, male and female, 
residing on the grid S = Z rf . Each site x G S contains two nests, one for 
the male particle and the other for the female particle. Each nest can be 
inhabited by at most one particle, either male or female. Let E = {0, 1} and 
F = E x E be the set of possible states at each site in S. For x £ S, we write 



where £ l (x) denotes the number (0 or 1) of male particles at site x and £ 2 (x) 
denotes the number of female particles at site x. We define the interaction 
neighborhood 

AT={0,yi,...,yjv} 

and the neighborhood of x 

M x = x + N. 

For example, N = {0, —1,1} if the interaction is nearest-neighbor on Z. Let 
Ci(x,m, £) denote the rate at which nest m (m = 1,2) of site x flips to state 
i (i = 0,1) and assume that Cj(x,m,£) depends only on the neighborhood 
M x , that is, 



while the birth rate ci(x,m,£) is positive only if both male and female 
particles can be found in J\f x . In this work, we take 5 = 1. Note that this 
simply means that the unit of time we take is the average lifetime of an 
individual. For example, the dioecious branching particle model which we 
will consider in Section 1.1 has 



ttx)=(e(x),f(x)) 



d (x, m, f ) = h i>m (£(x) ,€(x + yi),...,€(x + y N )) 
for some function hi :Tn : F N+1 The death rate cq is always constant, 






otherwise 




where 



n m/ (x,0 = \{zeM x :r'(x + z) = l}\ 
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that is, at rate A, each pair of male and female particles in M x give birth to 
a particle at nest m of site x if that nest is not already occupied. A more 
stringent condition, as in the particle model with rapid stirring which we 
will consider in Section 1.2, is to require both parent particles to reside at 
the same site, that is, 



Xn 1+2 (x,0, ifr(z)=0, 
0, otherwise, 



(3) ci(x,m,£) 
where 

ni +2 (s,0 =\{zeJ\f x :£ 1 (x + z) = 1 and f(x + z) = 1}|. 

This more stringent condition should not alter the behavior of the particle 
system if one allows a larger A than in (2), but it does help to simplify the 
analysis somewhat. 



1.1. Dioecious branching particle model. We first describe the model 
with birth and death rates as in (1) and (2), for which we will establish 
the existence of nontrivial stationary distribution(s) and, consequently, a 
phase transition (Section 2). The birth and death rates at site x depend 
only on £ in a neighborhood of x, therefore birth (and death) rates at x and 
y where N x C\N y = are independent. First, we restate the model: 

1. Birth. For each nest (x,m) and each pair (z±,Z2) G M x X N x such that 
q(z±) = 1 and ^ 2 {z2) = 1, where z\ and z 2 need not be distinct, with rate 
A, a child of (zi,z 2 ) is born into nest m of site x if (x,m) is not already 
occupied. 

2. Death. Each particle dies at rate 1. 

We can think of this particle system as a generalized spin system, gener- 
alized in the sense that the phase space at each site is {0, l} 2 rather than 
{0, 1}. One can refer to Chapter 3 of Liggett [5] for a detailed introduction to 
classic spin systems. We observe that the all-0 state [i.e., ^(x) =£ 2 (x) = 
for all x] is an absorbing state, therefore the probability measure that con- 
centrates only on the all-0 state is a trivial stationary distribution. We say 
that a stationary distribution is nontrivial if it does not concentrate only 
on the all-0 state. A major goal of this work is to establish the existence of 
nontrivial stationary distributions for various particle systems. 

The interacting particle system involving the birth and death mechanisms 
described above can be constructed using a countable number of Poisson 
processes, similar to the construction found in Chapter 2 of Durrett [2]. 
Define 

c* = sup Cj (x, m, £) . 
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We assume c* < oo. Let {T^' l ' m : n > 1} be the arrival times of independent 
rate c* Poisson processes and {U^' l,m :n > 1} be independent uniform ran- 
dom variables on [0,1]. At time t = T^' t,m , nest (x,m) flips to state i if 
jjx,i,m <- Ci (^ x ^ m ^ t _y c * anc i stays unchanged otherwise. 

Alternatively, one can explicitly write down the generators Q 1 and Q 2 
associated with the particle system with death rates (1) and birth rates (2) 
and (3), respectively, as follows: 



9 x m = E 

(x,rn)eSx {1,2} 



(4) 



r(x)(/(c-^, m )-/(0) 
+ E A£W(2)(i-rozO) 

y,zeAf x 



and 



(5) 



s 2 m) = E 

(x,m)eSx {1,2} 



/(0) 



x(/(e+^ 



f(0) 



where / is a function on F s (endowed with the product topology) with 
compact support (in this case, compact support implies that / depends on 
only finitely many sites) and S x m is a function on S x {1,2} that is one 
at (x,m) and zero elsewhere, and apply Theorem B3 in Liggett [6] to see 
that the closure of Q [in C(F S ) x C(F S ), with respect to the sup-norm 
topology] generates a Feller Markov process. Summarizing results from the 
two preceding paragraphs, we have the following theorem. 



Theorem 1.1. There exists a unique Feller Markov process £t con- 
structed as before with generator (4) or (5). 

The particle system £ with generator (4) or (5) is attractive in the sense 
that £ is monotonic in initial conditions. One can check that if £,o(x) < £o( 2; ) 
for all x £ S, where < denotes the partial order (0,0) < (0, 1), (1,0) < (1, 1), 
then £t( x ) — £,t(x) for all x and t. This is true since every birth or death 
event preserves <. For example, if Ct-(x) = (0,0), £,t-(x) = (0, 1) and at time 
t there is a male birth event at site x, then £t( x ) = (1)0) and £t( x ) = (L 1), 
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so the inequality £t( x ) < £t( x ) has been maintained. Similarly, one can check 
that the particle system £ is increasing in the birth rate A by coupling the 
random variables T^ l ' m and U^' l ' m involved in the constructions in the obvi- 
ous way. Because of this monotonicity, along with the existence of nontrivial 
stationary distributions for sufficiently large A and extinction for sufficiently 
small A which we will establish later in this work, we may conclude that 
there is a phase transition in the behavior of the particle system £. 

1.2. Description of the particle model with rapid stirring. If we add rapid 
stirring to the particle system, that is, we scale the integer grid 7L d by e and 
stir neighboring particles at rate s~ 2 in addition to performing the birth 
and death mechanisms, then the particle system converges to the solution 
of a reaction-diffusion PDE as e — > (see Theorems 8.1 and 8.2 in Durrett 
[2] and the beginning of Section 1.3 of this work). This PDE represents the 
mean-field behavior of the particle system and is usually easier to analyze 
than the particle system itself. As promised earlier, we will establish in Sec- 
tion 2 that there is a phase transition for the dioecious branching particle 
model (i.e., without rapid stirring), but obtaining any reasonable estimates 
on exactly where this transition occurs seems to be difficult. One advan- 
tage of adding rapid stirring mechanisms is that one can get a reasonably 
good idea of where the phase transition occurs in the rapidly stirred par- 
ticle model by analyzing the limiting PDE, or simulating this PDE on a 
computer. 

Moreover, this convergence establishes a connection between the particle 
model and PDE systems, which is of independent interest. Since many PDE's 
arise out of natural systems, this connection justifies the study of the PDE. 
The underlying stochastic system can also yield information about the PDE; 
for example, as we will see in Section 1.3, the monotonicity of the particle 
system will lead to the monotonicity of the PDE. Information about the PDE 
will similarly yield information about the particle model. In Section 4, we 
will establish a condition on the PDE that implies the existence of nontrivial 
stationary distributions for the particle system with sufficiently small e. 

For the particle models with rapid stirring, we work with S = eL d and 
denote the corresponding process by £ e . We also assume the birth rates 
in (3) and death rate 6 = 1, while the neighborhood N is nearest-neighbor: 
N = {u '■ \\y\\ = or e}. Here, we use the L 1 -norm: ||y|| = J2k=i \Vk\- in ad- 
dition to the transitions in the dioecious branching model, we introduce 
spatial movement of particles between neighboring sites called rapid stir- 
ring. We consider two rapid stirring mechanisms in this work, one called 
lily-pad stirring and the other called individual stirring: 

• Lily-pad stirring. For each x,y £ eL d with ||x — y||i = e, £ e (x) = (£, £,1 (x), 
£, £ ' 2 (x)) and £ e (y) = (£ e,1 (y), (, 6 ' 2 (y)) are exchanged at rate e~ 2 . 
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• Individual stirring. For each i G {1,2} and x,y G eZ^ with ||x — y\\\ = e, 
£ e '*(x) and £ £ ' ? (y) are exchanged at rate £~ 2 . 

Just as in the particle model without rapid stirring described in Sec- 
tion 1.1, one can construct the particle model with either lily-pad stirring or 
individual stirring using a countable number of Poisson processes. Alterna- 
tively, one can write down the generator explicitly and again apply Theorem 
B3 in Liggett [6] to establish: 

Theorem 1.2. Let S = eL d and M = {y: \\y\\ = or e}. There exists a 
unique Feller process £t with generator Q L for the particle model with lily-pad 
stirring or generator Q 1 for the particle model with individual stirring, 

(6) G L f(Z) = G 2 f(0+ E e" 2 (/(^)-/(0), 

x,y£S,x£Af y 

(?) g I f(o=Q 2 f(0+ E e- 2 v(n {x ' m) " {y ' m) )-fm, 

me{l,2},x,yeS,x£Af y 

where 

£(z,m'), ifz^x,y, 

£(2,771'), ifz = y, 

£(y,m'), ifz = x 

£(z,m'), if(z,m')^(x,m),(y,m), 
£(x,m), if (z,m') = (y,m), 
€(y,m), if (z,mf) = (x,m). 

For lily-pad stirring, instead of thinking of a site that consists of two 
nests, as in the dioecious branching model, we can view each site as having 
for states in 

F = {0,1} 2 = {(0,0),(0,1),(1,0),(1,1)}. 

We restate the dynamics of the particle model in terms of these four states. 
At any site x G eL d ', only the following transitions are possible: (0, 0) <-> (0, 1), 
(0,1) (1,1), (0,0) <-> (1,0) and (1,0) <-> (1,1), that is, only one particle is 
born or dies at a particular time. The rates of these transitions are as follows: 

c m (x,n = h iff (2) = (0,1) or e (x) = (1,0), 

c { o,i)(x, e ) = c m ( x , e) = 1, if r (x) = a, i), 

C(o,i)(x,£ e ) = C(i,o)(x,£ £ ) = \ni +2 (x,£ £ ), 

if f(2) = (0,0), 

c { i,i)(x,e ) = An 1+2 (x,f ), if e (x) = (0, 1) or f (x) = (1,0). 
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The difference between these two stirring mechanisms is that lily-pad 
stirring forces male and female particles at a site to move together, but in- 
dividual stirring allows independent movement of male and female particles. 
Every exchange of particles, in both lily-pad stirring and individual stir- 
ring, is monotonicity preserving, thus neither stirring mechanism disrupts 
the monotonicity property of the particle system. 

1.3. Convergence to a PDE for lily-pad stirring. Consider the particle 
system with lily-pad stirring and its generator given by (6). For i 6 F, if we 
define 

ui(t,x)=p(e t (x)=i), 

then Theorem 8.1 in Durrett [2] (or its generalization, Theorem 3.1 of this 
work) shows that if g%{x) is C 1 and uf(0,x) = gi(x) for all i, then 

(8) Ui(t,x) = limufft, x) 

£—♦0 

exists and satisfies the following system of PDEs: 



(9) 



dt 

du(o,i) 
dt 

d«(i,o) 
dt 



& u (o,o) + «(o,i) + n (i,o) - 4^(0,0)^(1,1)* 
Au (0 ,i) + - n (0 ,i) + 2\d(u {m - u (0j i) 

Au ( i,o) + «(i,i) - "(i,o) + 2 M w (o,o) - W(i,o))n (1 , 1 ), 

Au (l,l) ~ 2u (l,l) + 2 M n (0,l) + «(1,0) )«(!,!)■ 



dt 

In order for the limit in (8) to make sense, we extend the definition of 
uf(t,x) to all x G M. d by requiring uf(t,x) = uf(t,e\_x/s\) where, [q\ denotes 
the integer part of q. Obviously, must lie in [0, 1] for all i and t since 
it is a limit of probabilities. The reaction terms of the equation involving 
du(p$)/dt say roughly that sites enter state (0,0) at rate W(o,i) +"(i,o)> that 
is, when the only male or female particle at the site dies; sites leave state 
(0,0) at rate 4Adti( ,o)' u (i,i)) that is, when a particular site has no particles 
and a particle (male or female) is born by a pair of male and female particles 
at one of the 2d neighboring sites. 

We also observe that, strictly speaking, one should require the initial con- 
dition of (9) to be C 2 , for otherwise, the PDE system may not make sense 
at t = 0. This issue can be remedied by considering the corresponding inte- 
gral equation, as in equation (15.1.2) of Taylor [7]. Furthermore, Proposition 
15.1.2 of Taylor [7] shows that solutions to (9) are C°° at any t > as long 
as the initial condition is C 1 . Thus we only require our initial conditions to 
be C 1 from this point on. 
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We transform the parameter space of the 3-dimensional system in (9) to 
obtain a monotone 2-dimensional system, which will be easier to analyze. 
First, define v = U(o,o)> v i = u (o,i) + n (i,o) and v 2 = «(i,i), (vo,vx,v 2 ) then 
satisfies 

— = Av + v i - 4Xdvov 2 , 
ot 

= Avi + 2v 2 - vi + 2Xd(2v - vi)v 2 , 

dv 

-p— = Av 2 - 2v 2 + 2\dviv 2 , 
Ot 

where vq + vi + v 2 = 1. The above system can be written as the limiting 
PDE under rapid stirring of another particle system £ £ , still on S = eZ d , 
with state space F = {0, 1, 2} and transitions <-> 1 and 1 2 at rates 

Co (x,c £ ) = i, if c £ (^) = i, 

Cl (x,C e ) = 2, if C £ (^) = 2, 

Cl (x,C e ) = 2An 2 (x,C e ), if C e (^) = 0, 
c 2 (x,C £ ) = \n 2 (x,C), if( e (x) = l, 



where 



n 2 (x,e) = \{zeN:C e (x + z)=2}\. 



One can interpret C £ ( x ) as the total number of ^-particles at x and check 
that monotonicity in initial condition holds for £ e . Let (t>f (0, x), v 2 (0, x)) = 
(gi(x), g 2 (x)) and (vf (0,x),v 2 (0,x)) = (gi(x),g 2 (x)) be two sets of initial 
distributions such that g 2 < g 2 and g\+g 2 <gi+g 2 , with g 1 , g 2 , g 1 + g 2 , g x , 
52; gi + 52 all lying in [0, 1]. It is then possible to set up two initial conditions 
Co and Co such that P((q(x) =i) =vf(Q,x) and P(Q(x) = i) =vf(0,x), i = 
1,2, and Co( x ) — Co( x ) holds for all x and u. The monotonicity property of 
C e implies that Ct( x ) — Ct( x ) f° r all t and x and therefore 

P(C!(x) > 1) < P(Cf (x) > 1) and P(C t £ (x) > 2) < P(£(s) > 2), 

that is, for all t and x, 

v 2 (t,x) < v £ 2 {t,x), 

vl(t,x) + v 2 (t,x) < vf(t,x) + v 2 (t,x). 

We now transform the parameter space a second time by defining (u,v) = 
(v\ +v 2 ,v 2 ) and writing c= 2\d. (ut,Vt) is then monotone in initial condi- 
tions since vf(t,x) — ► Vi(t,x). We observe that u can be interpreted as the 
density of occupied sites (where either one or both nests are occupied) and 
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v as the density of doubly occupied sites (where both nests are occupied). 
Straightforward calculation shows that (u, v) satisfies the following system: 

du 

— = An + (2c(l - it) + l)v - u, 
(10) * 

— = Av + (c(u - w) - 2)u. 
at 

We summarize this paragraph in the following lemma. 

Lemma 1.3. The PDE system (10) is monotone in initial conditions 
that lie in 

K = {{u,v):0<v<u<l}, 

that is, if there are two initial conditions (uo,vo) £ 1Z and (uq,vo) E 1Z, with 
uq < uq and vq < vo everywhere, then u t < u t and v t < v t everywhere, for 
all t. Furthermore, both (ut,vt) and (ut,vt) lie in 1Z for all t. 

In Section 4, we will analyze (10) to establish the following result. 

Theorem 1.4. If X is sufficiently large and e is sufficiently small, then 
there exists a nontrivial translation invariant stationary distribution for the 
dioecious branching particle model with lily-pad stirring with generator (6). 

1.4. Convergence to a PDE for individual stirring. We consider the par- 
ticle system with individual stirring and with its generator given by (7). 
Unlike lily-pad stirring, Theorem 8.1 in Durrett [2] cannot be directly ap- 
plied to obtain convergence to a PDE system for individual stirring. We can, 
however, follow the ideas used in the proof of that theorem to establish a 
corresponding result, Theorem 3.1. For the process £ £ with generator (7) 
and i £ E, define 

u{ m {t,x) = P{£ £ t {x,m) = i). 

Theorem 3.1 then implies that if gi jTn :W^ [0,1] is C 1 and uf m (Q,x) = 
9i,m( x )i then Ui^ m (t,x) = lim e ^o u \ m (t, x) exists and satisfies the following 
system of PDEs: 



dt 
duio 



- = Aui ; i - ui 5 i + 2\d(l - ni, 1)^1,1^1,2, 
Aui,2 - ui, 2 + 2Xd(l - 1/1,2)^1, 11*1,2, 



dt 

where uq 1 + u\ t i = UQ t 2 + = 1- Notice that if we start with a symmetric 
initial condition, that is, g.;,i = ^2, then the solution to the above PDE is 
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also symmetric. And if we define u = ux t i = Ui t 2, then we obtain the following 
PDE for u: 

(11) — = At* + /(«), f{u) = -u + 2\d(l-u)u 2 . 

This PDE has been analyzed in Durrett and Neuhauser [3] as their sexual 
reproduction model (Example 3 on page 291). In fact, it is not difficult to 
see that if U\ \ = U\2, then choosing the "father" from the male population 
is exactly the same as choosing the "father" from the female population, 
hence it is quite natural for this reduction to occur. Theorem 4 of Durrett 
and Neuhauser [3] states that if 2 Ac? > 4.5 and e is sufficiently small, then 
their sexual reproduction model has nontrivial stationary distribution(s). 
Although this theorem does not directly apply to our particle system £ e 
with two types of particles because of the difference in stirring mechanisms, 
one can nevertheless work through the proof of Lemma 3.3 of Durrett and 
Neuhauser [3], while making obvious changes, to establish the following, 
similar, result. 

• Let < p\ < po < 1 be the two nonzero roots of f(u). Define (3 = (po — 
pi)/10 and Ik = 2Lke\ + [— L,L) d . If e is small, L is large, and £ e (0) has 
density at least p\ + (3 of both male particles and female particles in Iq , 
then for sufficiently large T, with high probability, £ e (T) will have density 
at least po — (3 in I\ and 

This result can then be fed into a comparison argument, comparing the 
particle system with oriented percolation, as on page 312 of Durrett and 
Neuhauser [3] or in the proof of Theorem 4.3 in Durrett [2], to establish the 
existence of nontrivial stationary distribution(s) for the particle system £ £ 
under individual stirring with sufficiently small e. We then have the following 
theorem. 

Theorem 1.5. If X is sufficiently large and e is sufficiently small, then 
there exists a nontrivial translation invariant stationary distribution for the 
dioecious branching particle model with generator (7). 

1.5. Discussion. We can numerically solve PDEs using standard finite 
difference methods to obtain rough values of A c for d = 2, where A c is, as 
in Durrett and Neuhauser [3], the critical value of A necessary for longterm 
survival in the rapid stirring limit (e — > 0). We recall that (10) and (11) 
are the rapid stirring limits of the particle systems with generators Q L 
and G 1 in (6) and (7), respectively, both with birth-death mechanisms 
described by Q 2 that require parents to be at the same site. We can re- 
place Q 2 with Q 1 (which only requires one parent to be in the neighbor- 
hood of the other parent) in (4) and obtain slightly different generators 
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Individual stirring 


Lily-pad stirring 


G 1 


A c = 0.225 


A c G [0.271,0.272] 


9 2 


A c = 1.125 


A c € [1.114, 1.115] 



G L f(0 = S l f{i) + Ex,yeS,xeM y e- 2 (f(^ y ) ~ HO) and G 1 f(0 = 0^(0 + 
E me{ ii}, x , ye s^N y Z~Hf(&> m) ~ (y > m) ) ~ HO)- The rapid stirring limit of 
Q 1 is exactly the PDE in (11) but with 2d replaced by 2d(2d + 1) because 
there are 2d(2dA- 1) "potential parent-nest pairs" involved in Q l . The rapid 
stirring limit of Q L satisfies the following PDE: 

ft = A «(o,o) - 2A(2(i) 2 u ( o,o)(w(o,i) + «(i,i))(«(i,o) + u (i,i))> 

dt = Au( 0j i) - A(2d)(2d + l)u( ,i)(u( ,i) + «(i,i))(u(i,o) + «(i,i)) 
+ X(2d) 2 u (u {0A) + U(i,i))(u(i j0 ) + «(i,i)), 

dun o) 

( 12 ) g t = A«(i i0 ) - X(2d)(2d + l)«(i |0 )(«(o,i) + «(i,i))(«(i,o) + «(i,i)) 
+ A(2d) 2 u (M (0 ,i) + n (lj i))(u (li0) + 1^(1,1)), 

Au(i,i) - X(2d)(2d + l)(u( 0) i) + «(i,o)) 

X ( w (0,l) + «(i,i))( u (i,o) + u (i,i))- 

For individual stirring, the limiting PDE is (11) and Theorem 4 from 
Durrett and Neuhauser [3] shows that A c = 4.5/2d = 1.125 for C/ 7 and A c = 
4.5/2d(2d+ 1) = 0.225 for For lily-pad stirring, we do not have a readily 
available theorem to tell us the exact value of A c . We obtained the range of 
values of A c in the rightmost column of Table 1 by simulating (10) and (12) 
using progressively finer grids and stopping once changes in the estimates of 
A c become smaller than 10 -3 . 

Curiously, free movement of all individuals (individual stirring) seems 
to work better (to ensure survival at least) when mating occurs between 
all pairs of individuals in the neighborhood, as in Q 1 , but restricting male 
and female individuals to move together (lily-pad stirring) seems to work 
better when mating occurs only between male and female individuals at the 
same site, as in Q 2 . A possible explanation is that matching of movement of 
individuals and the mating strategy helps survival: when the mating strategy 
is to mate with any individual of the opposite sex in the neighborhood, 



d"(i,i) 
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free movement of all individuals helps everyone to find mating partners 
more easily; but when the mating strategy is to mate with only individuals 
at the same site, free movement of all individuals only helps to break up 
"marriages." 

In the single-sex scenario where the birth and death mechanisms are the 
same as those of the contact process with rate of infection A and rate of 
recovery 1, we obtain the following PDE in the rapid stirring limit: 

du 

— = Au + f(u), f(u) = -u + 2Xd(l-u)u, 

which has A c = l/2d (i.e., A c = 0.25) by Theorem 2 from Durrett and 
Neuhauser [3] . This A c is roughly the same as the A c for birth-death mecha- 
nisms Q l , although much smaller than the A c for birth-death mechanisms Q 2 
since the birth mechanism in Q 2 is much more restrictive. With birth-death 
mechanisms Q 1 and d = 2, the density of individuals (male or female) at 
equilibrium when A is only slightly higher than the A c is roughly 0.66 for in- 
dividual stirring and 0.79 for lily-pad stirring. Thus, in the two-sex scenario 
Q\ although each individual can only give birth if there is any individual of 
the opposite sex in the neighborhood, finding a mate in the neighborhood 
should not prove to be a problem and it is not terribly surprising that the 
A c for Q 1 with individual stirring is smaller than the A c for the single-sex 
scenario. 

In the remainder of this paper, we will establish various results as promised 
in this section. In Section 2, we prove a few results on the dioecious branch- 
ing particle system without rapid stirring. In Section 3, we present, without 
proof, a convergence theorem for particle system with individual stirring. Fi- 
nally, in Section 4, we prove Theorem 1.4 by establishing a condition on (10), 
similar to condition (*) on page 180 of Durrett [2]. 

2. Results on the dioecious branching particle model. In this section, 
we assume the model with generator (4) where we take S = that is, 
the particle system with birth and death mechanisms, but no stirring. We 
briefly restate the model: the rate at which nest m of site x flips to state i, 
Ci(x,m,g), is 

1, if£ m (x) = l, 
0, otherwise, 



Cl(x ' m ' e) = \0, otherwise, 

where n m /(x,£) = \{z £ N x : £, m ' (x + z) = 1}| and M x contains the site x and 
its 2d nearest neighbors. The goal is to establish the existence of a phase 
transition. 
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2.1. Existence of stationary distributions. We first establish that station- 
ary distributions exist. This is a generalization of Theorem 2.7 in Durrett 
[2] or Theorem III. 2. 3 in Liggett [5]. We use a method along the lines of 
Theorem 2.7 in Durrett [2]. Define 

£o(%) = (1) 1) f° r an x 

and let & be the process started with initial condition £ - Let T f /(£o) = 
E f(£t) be the semigroup corresponding to the particle system. Tt is then 
a Feller semigroup by Theorem 1.1. We begin with a lemma. 

Lemma 2.1. For any A,E> C 5 = 7L d , the junction 
(13) t»P(£(x) = VsG A,g(y) = Vy G B) 

is increasing. 

Proof. Let ao = £s an d Aj = £s f° r an arbitrary fixed s. Then £o( x ) — 
cto(x) and £o(x) > @o(x). Let {at, (it) be the state at time i of the particle 
system that started with initial condition (ao>A))- Then, by the fact that 
the particle system is monotone in initial conditions, we have 

$(x)>a t (x) and $(x)>p t (x) 

for all t and x. The Markov property of £ then implies that the function 
in (13) is increasing in t. □ 

Theorem 2.2. As t—* oo, & converges weakly to t;^. The limit is a 
stationary distribution that stochastically dominates all other stationary dis- 
tributions and is called the upper invariant measure. 

Proof. Let A and B be arbitrary subsets of S. For C = {x±, . . . , x m } C 
S and D = {yi, . . . , y n } C S, we write 

P(£(z) = Vz G A£?H =0Vw£ B.g (x) = 1 Vx G C,^(y) = 1 Vy G £>) 

(m+n \ 
u E *y 

where 

^ = ( z ) = Vz G A U {xj, g(t») = 0VweB} if i = 1, . . . , m 

and 

#i = U7 (*) = Vz G A, £? H = Vw; G B U {2/i_ m }} 

if ? = m + 1, . . . , m + n. 
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We can use the inclusion-exclusion formula on PiXSj^y 1 Ei), that is, 

/m+n \ m+n 

p ( U E i ) = E p (^)-E p (^ n ^) 

\ i=l I i=l i<j 

+ - + (-i) m + n+1 P(E 1 n — nE m + n ). 

Every term in the above expansion is in the form 

P$(z) = Vz G ;g(w) = Vw G ••), 
which is increasing in t by Lemma 2.1. Therefore, 

P(g(z) = Vz G A£?H = Viu G fl,g(a;) = 1 Vx G C,g(y) = 1 Vy G D) 

converges for all A, B, C and D where C and D are finite, hence all finite- 
dimensional distributions converge. Thus, a weak limit (£^,£^0) exists and 
since T t is a Feller semigroup, (£^'£00) ^ s a stationary distribution. We can 
also easily see, via a monotonicity argument, that {Clc^oo) dominates all 
other stationary distributions. See the proof of Theorem 2.7 in Durrett [2] 
for a simpler version of this type of argument. □ 

2.2. Extinction for sufficiently small A. 

Theorem 2.3. If \\J\f\ 2 < 1, then the particle system £ with genera- 
tor (4) has no nontrivial stationary distribution. 

Proof. We compare a modification of the particle system £ with the 
contact process. Theorem 2.6 of Durrett [2] states that if a\N\ < 1, where a 
is the rate of infection, then the contact process has no nontrivial stationary 
distribution. 

We modify the birth rates ci(x,m, £) in (2) and (3) of the particle model 
to obtain 

(A|.Af|ni(x,£), if£ 1 (^) = 0, 
c' 1 (x,m,0 = < X\Af\n 2 {x,0, if£ 2 (x) = 0, 
1 0, otherwise, 

that is, birth of male (female) offspring no longer requires the presence of 
female (male) parents in the neighborhood. We denote this modified process 
£ = (£ 1 ,£ 2 ). The result of the modification is that £ x and £ 2 are now decou- 
pled and £ J behaves in the same way as the contact process with birth rate 
a = \\J\f\. Furthermore, the modified process (£ 1 ,£ 2 ) can also be constructed 
using the Poisson processes {T^' vm :n > 1} and {U^' l ' m :n > 1}, just as in 
Section 1.1. Since c^(sc,m,£) > ci(x,m,£) for all x, m and £ [with the def- 
inition in either (2) or (3)], (£ x ,£ 2 ) stochastically dominates the original 
process (£ 1 ,£ 2 ). If a\M\ < 1, then £ has no nontrivial stationary distribution 
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and converges weakly to the all-0 state as t — > oo for any initial con- 

dition. Thus, (£t ,£t ) also converges to the all-0 state for any initial condition 
if a\N\ = \\M\ 2 < 1, as required. □ 

2.3. Survival for sufficiently large A. We use the idea of Chapter 4 of 
Durrett [2], that is, we compare the particle system to an oriented perco- 
lation process, see Durrett [2] for more details on the oriented percolation 
process. A particularly useful result that will be used in the proof below is 
Theorem 4.2 of Durrett [2]. 

Theorem 2.4. Let Wfi be an M-dependent oriented percolation process 
with density at least 1 — 7, starting from the initial configuration Wq in 
which the events {x £ Wq }, x £ 2Z, are independent and have probability p. 
Ifp>0 and 7 < 6" 4 ( 2M+1 ) 2 , then 

limjrfP(0€W&)> jjj. 

This theorem shows that if the density of open sites 1 — 7 is sufficiently 
close to 1 and we start with a Bernoulli initial condition for Wq, then the 
probability that is wet at time t does not go to as t — > 00, hence the 
upper invariant measure is nontrivial. 

Theorem 2.5. If X is sufficiently large, then the particle system £ with 
generator (4) has a nontrivial stationary distribution. 

Proof. We follow the method of proof as in Chapter 4 of Durrett [2]. 
First, we describe a construction of the particle system £ that is more spe- 
cialized than the one given in Section 1.1. Recall that S = "L d and N = 

{x E S:\xx\ H h \x d \ = or 1}. Let m £ {1,2}, x,y,z£ S, {R% m ,n> 1} 

be independent Poisson processes with rate 1 and {T^' m ' y ' z , n > 1}, with 
y,z gM x , be independent Poisson processes with rate A. At time R^ m , any 
particle residing at (x,m) is killed. And at time T^' m ' y ' z , a particle is born 
at (x,m) if (x,m) is not already occupied and nests (y, 1) and (z, 2) are both 
occupied. 

We will select an event G^ , measurable with respect to the filtration 
generated by all of the Poisson arrivals {i?^' m } and {T^' m ' v ' z } at 

zE{(-l,0,...,0),(0,0,...,0),(l,0,...,0)} 

during the time interval [0,t). For any 7 > 0, regardless of how small, there 
exists A and T and an event G^ with 



P(G ?0 )>l-7, 
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so that on G 5o , if £ (0, 0, . . . , 0) = (1, 1), then f T (l, 0, . . . ,0) = £t(-1, 

0. . . . , 0) = (1, 1). One can achieve this by choosing T so small that the prob- 
ability of any death occurring at any nests at sites (—1, 0, . . . , 0), (0, 0, . . . , 0) 
and (1,0,..., 0) is less than 7/2; one can then choose A sufficiently large 
so that the probability of having birth events from ((0, 0, . . . , 0), 1) and 
((0, 0, . . . , 0), 2) to each of the four nests at sites (—1, 0, . . . , 0) and (1,0,..., 0) 
during [0, T) is larger than 1 — 7/2. In other words, if we define the event 

G^ = {There are no death events during [0, T) 

at sites (-1, 0, . . . , 0), (0,0,..., 0) 

or (1,0,..., 0); and there are birth events from ((0, 0, . . . , 0), 1) 

and ((0, 0, . . . , 0), 2) to each of the four nests at sites (— 1, 0, . . . , 0) 

and (1,0,..., 0) during [0,T)}, 

then Cr£ satisfies the requirement and P{G^ Q ) > 1 — 7 for some A and T. 
G^ is the "good event" that will ensure male and female particles get born 
at sites x — 1 and x + 1 provided site x is inhabited by both a male and a 
female particle. 

Using this "good event" Gg , we can construct an oriented percolation 
process of density at least 1 — 7 that is stochastically dominated by the 
particle system £, such that existence of nontrivial stationary distribution(s) 
for the oriented percolation process, as provided by Theorem 2.4, implies 
existence of nontrivial stationary distribution for the particle system. This 
part of the proof follows the proof of Theorem 4.3 in Durrett [2]. We therefore 
omit the details and instead refer interested readers to the proof of Theorem 
3.3.2 in Yu [9]. □ 

3. Convergence theorem for individual stirring. In this section, we present 
the convergence result for the individual stirring model, as promised in Sec- 
tion 1.4. We work in a slightly more general setting and consider random 
processes 

Ct:eZ d x {1,2,...,M}^{0,1,..., K -1}. 

We call each x € eL d a site and each (x,m) G eTL d x {1,2, . . . ,M} a nest. 
There are M nests at each site. We think of the set of spatial locations 
Z d x {1, 2, . . . , M} as consisting of M floors of Z d . Let 

J\f = {Q,ey 1 ,...,ey N } 

be the interaction neighborhood of site 0. The process £| evolves as follows. 

1. Birth and death. The state of nest (x, m) flips to i, i = 0, . . . , k — 1, at rate 

Ci(x, m, £) = h i>m (£(x, m),£(x + ezi, mi), . . . , £(x + ez L ,m L )), 
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where L is a positive integer, z±, . . . , zl £ M , mi, . . . ,niL € {1, 2, . . . , M} 
and 

/i Jim :{0,l,...,K-l} L+1 ^]R + 

with hi, m (i, . . .) = 0. 
2. Rapid stirring. For each m £ {1, 2, . . . , M} and x, y £ eZ d with ||x — y||i = 
e, £ e (x,m) and £, e (y,m) are exchanged at rate e -2 . 

This individual stirring model differs from the lily-pad stirring model de- 
scribed in Section 1.2 in that the stirring actions between corresponding 
nests at neighboring sites are now independent. More specifically, exchanges 
are allowed between neighboring nests on the same floor only, that is, be- 
tween (0,1) and (e,l), but not between (0,1) and (e, 2). 

As an example, for d = l, in the particle model with individual stirring 
with generator (7), we have k = 2, M = 2, L = 4, N = {0, — e,e}, 

1, if £(x, m) = 1 
0, otherwise, 



A(£(x - e, l)£(x - e, 2) + £(x + e, 2)£(x + e, 1)), 
ci(x,m,£) = { if £(x,m) = 0, 

0, otherwise. 

In particular, we should define 

(zi,mi) = (-1,1), (22, m 2 ) = (-1,2), 

(z3,rn 3 ) = (1,1), (z 4 ,m 4 ) = (1,2) 

and /ij = as 

^o("0, "1,02, 03,04) = "0, 

/ii(ao,«i, «2, «3, 04) = A(ai«2 -I- 0304) (1 — qo). 

Theorem 3.1. Suppose {^(x,m),(x,m) e eZ d x {1,2, ... ,M}} are in- 
dependent and let u £ im {t,x) = P(£f(x,m) = i). If uf m (0,x) = g^ m {x) and 

gi^ m :M d — ► [0, 1] is C 1 for all i and m with J2i9i,m = 1» then, for any smooth 
function (f> with compact support, as e — > 0, 

(14) e d 0(y)l{ef(y,m)=i} ~> 4>{y)ui, m (t,y)dy in probability, 

y eeZ d 

where Ui im (t, x) zs i/ie bounded solution of 
dui , 



^^i,m "I - /i,m(^), ^i,m(0, x) — £/j jm (x), 



at 

/i,m(«) = (ci(0, m, e)l(e(0, m) / »)>« " E( C J (°' m ' 01(^(0, m) = i)) u 
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and {(/)(£)) u denotes the expected value of <fi(£) under the product measure in 
which state j at nest m has density Uj >m , that is, ^{x,rn), with x G eWj and 
1 < m < M , are independent, with P(£(x,m) = j) = Uj >m . 

We refer interested readers to Theorem 4.0.5 of Yu [9] for its proof, which 
follows the approach used in the proof of Theorem 8.1 in Durrett [2]: first, 
a dual process is defined for the particle system, then this dual process is 
shown to be almost a branching random walk that converges to a branching 
Brownian motion as e — ► 0; furthermore, two different duals are asymptoti- 
cally independent of each other. This asymptotic independence implies (14), 
as well as the convergence of the particle system itself to Ui >m . 

In the proof of Theorem 3.1, the definition of the dual is the only part 
that differs from the proof of Theorem 8.1 in Durrett [2] in any significant 
way: each nest at a site requires a dual process and it is affected by any birth 
or death events that happen to any nests at that site, but dual processes for 
any two nests (even at the same site) are asymptotically independent. 

4. Invariant stationary distribution for lily-pad stirring. In this section, 
we establish the existence of nontrivial stationary distribution(s) of the par- 
ticle system with lily-pad stirring (as promised by Theorem 1.4) by showing 
that for sufficiently large c, the solution to 

du 

— = Au + (2c(l -u) + l)v - u, 

(15) 61 
dv 

— = Av + (c(u-v)-2)v, 

with initial condition uq = /, vq = g, f > g satisfies the following condition: 

Condition (*). There are constants < D\ < d\ < d 2 < D 2 < 1, M and 
T such that if v (x) G (D 1 ,D 2 ) for x G [-M,M], then v T {x) G (d 1 ,d 2 ) for 
x£ [SM,3M]. 

According to Chapter 9 of Durrett [2], this is a sufficient condition for 
the existence of nontrivial invariant stationary distribution(s) for the par- 
ticle system with sufficiently fast stirring, so Theorem 1.4 will follow once 
Condition (*) is established. Recall that Theorem 2.5 establishes that the 
dioecious particle model without rapid stirring has a nontrivial stationary 
distribution if the birth rate A is sufficiently large. If one works through 
the proof, however, one will find that "sufficiently large" in that argument 
means that A is larger than a number on the order of 6 100 , which is not too 
informative as to where exactly the critical A for the phase transition is. On 
the other hand, in the model with rapid stirring, one can obtain a far better 
idea of the range of A for which Condition (*) holds. 
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In this work, we also establish Condition (*) for sufficiently large c (recall 
that c = \d), but here, "sufficiently large" means that c is only larger than 
a number on the order of 10 4 . We assume dimension d = 1; extension to 
d > 1 is straightforward. The proof consists of two parts: the first part, 
Section 4.1, establishes the existence of constants d\ and D\ and the second 
part, Section 4.2, establishes the existence of constants d<i and D2; the second 
part will be easy once the first part has been established. 

4.1. Lower bounds: Existence of d\ and D\ in Condition (*). Let 

K = {(u,v):0<v<u< 1} 

and 77 be a vector field on TZ that generates a flow (J-^)t>o on TZ, that is, if 
(uq,vq) £ TZ, then J^(uq,vo) £ TZ for all t > 0. We assume T 1 ^ is monotone 
(in initial conditions), that is, it preserves the partial order on TZ given by 

(ui,vi) < (ii2,V2) <==^ (tii < v-2 and v\ < V2)- 

In other words, if {u\,V\) < (^2,^2), then Fi(u\,v\) < .^(t^, ^2)- The sce- 
nario for 77 that we consider is the following: the ODE system 

du . dv , . 

— = m (u,v), — = m {u,v) 

has a stable fixed point -P+ close to the top corner of TZ with a relatively 
large basin of attraction, but (0,0) is also a stable fixed point (with a much 
smaller basin of attraction), and there is another unstable fixed point P_ 
lying "between" these two stable ones; see Figure 1 for two examples of the 
vector field 77 found in (15). 

Theorem 9.2 in Durrett [2] establishes Condition (*) for a specific predator- 
prey system with phase space {0, 1,2} at each site. The critical fact used in 
the proof is that the associated ODE system has only one interior equilib- 
rium point and has a global Lyapunov function. The phase portrait of the 
ODE associated with (15), however, shows that it has two interior equilib- 
rium points, P+ and P-, where P- is always a saddle point. Hence, we have 
a scenario where there is little hope of finding a global Lyapunov function 
and, in general, even finding an explicit Lyapunov function that works inside 
the domain of attraction of P+ is difficult. 

For general scalar reaction-diffusion equations (i.e., where the reaction 
terms are 1-dimensional), Chapter 15.4 of Taylor [7] provides an overview 
of methods and results. Convergence results in the scalar case, such as Con- 
dition (*) in Chapter 9 of Durrett [2], can be established using techniques 
found in Fife and McLeod [4], just as in Durrett and Neuhauser [3]. For mul- 
tidimensional reaction-diffusion equations, results regarding the existence of 
traveling wave solutions are more limited. Theorems 1.1, 1.2 and 4.2 of 
Chapter 3 of Volpert, Volpert and Volpert [8] provide existence results for 
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traveling wave solutions for certain classes of monotone reaction-diffusion 
systems. Indeed, Theorem 1.1 of Chapter 3 of Volpert, Volpert and Volpert 
[8] applies to equation (13), but estimating the speed of the wave [which 
is essential for ensuring that Condition (*) expands rather than shrinks] 
is still nontrivial and must be done on a case-by-case basis. An alterna- 
tive approach to establishing Condition (*) for equation (13) may be to use 
Theorem 1.1 of Volpert, Volpert and Volpert [8], which implies existence of 
traveling waves for (13), and then to try to obtain estimates for the speed 
of the wave (probably also a result involving sufficiently large c). With this 
estimate on the wave speed, one may then be able to use generalization of 
techniques in Fife and McLeod [4] to establish a convergence result. 

The method we use to establish Condition (*) for solutions of (15) is much 
more elementary and only requires the monotonicity property. Thus, it even 
applies to cases where the existence of traveling wave solutions is not known. 

Consider the PDE system in one spatial dimension, 



(16) 



du 
dt 



Au + rji(u,v), 



dv 
&t 



Av + rj2{u, v). 



We first define the shape of the initial conditions (no, vq), which is a smoothed 
indicator function of a suitable interval. Let 



(17) fo(x) 




if x G [-L + l,L- I], 

if xe [-L-l,-L + l], 

ifx£ [L — l,L + l] } 

if x e (oo, -L-l]U[L + I, oo) 




(a) c = 5 (b) c = 25 

Fig. 1. Phase space of the ODE: 771 = on 71 and 772 = on 72. 
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where fo S H 



(18) 



and 



h(x) 



x + l 



I 



x < —I, 
-l<x<0, 

0<x<l, 
x > L 



see Figure 2. In the above definition, the choice of L is arbitrary, provided 
L > I, but later, we will choose I small such that |A/o| is large in [— L — 
l,-L + l]U[L-l,L + l]. We call the intervals [-L-l,-L + l] and [L-l,L + l] 
the transition regions. We observe that both h and b! are continuous at x = 0, 
with 



/i"(x) 



1 

T 2 ' 



1 



if -Z < x < 0, 



if < x < L 



(19) 



|A/ | < 



so the graph of h in the plane is symmetric about the point (0, 1/2) and also 

1 
P 

everywhere. 

We assume the solution (u(t),v(t)) to (16) starts with initial condition 
(do/oj^o/o)- We would like to show that the interval over which the v- 
coordinate of the solution to (16) is > bo expands with time. Thus, intu- 
itively, we would like to see the ^-coordinate of (u(t),v(t)) increase, at least 
when v(t) is larger than some threshold but smaller than ir v (P + ), where 
7r v (u' ,v') =v'. If this were the case, we could use the v -coordinate as a "par- 
tial Lyapunov function" inside a subset of the basin of attraction of P+. 
Unfortunately, this does not always hold, as is the case considered in this 
work where rj = ((2c(l — u) + l)v — u, (c(u — v) — 2)v) and r]2(u,v) < when 
u — v is close to 0. 

To overcome this difficulty, we define a convex family of nested subsets 
of 1Z into which the flow Tn contracts at a sufficiently large rate. The up- 
per boundary of these nested subsets consists of lines u = v and u = 1, 
while the lower boundary consists of curves from the family {"fe}e£[o,e ] ) 
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where 70 = (1 + 0)7 n 1Z and 7 is a nonincreasing (i.e., with nonincreas- 
ing ^-coordinate when parameterized according to increasing it-coordinate) 
curve whose exact definition may change depending on the exact form of 
77. Assume, furthermore, that there is a family of continuous mappings 
<Pq : [0, 6>o] x 70 — > 70 such that T n satisfies the following. 

Assumption 4.1. There exist e, K\ , AT 2 > (K\ is large and is small) 
and sufficiently small so such that for 6 £ [0, 6>o], (oo> &o) £ 70 and s £ [0, so], 
we have 



where (a s ,b s ) = tp e (s,(a ,b )). 

For each s, the mapping ipg(s,-) maps points on the curve 70 to other 
points still on 70 . Under T r ^ the upper part (the part above the horizon- 
tal line v = e) of the line segment {(aao,abo) :0 < a < 1} is pulled above 
(aa s ,afr s ), while its lower part is not pushed too far beneath (aa s ,ab s ); see 
Figure 3 for a schematic illustration. Note that ((1 + Kis)a s , (1 + K\s)b s ) 
again lies on a curve 7^, where 6 > 9 + 5s and 5 depends on K\, K2 and the 
geometry of 7. 




if afro > e, 
if afro < £) 



1 




M <^(«o,&o) 



f 







u 



1 



Fig. 3. Illustration of Assumption 4-1- 
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Proposition 4.2. Assume that J-„ satisfies Assumption 4.1 with 
min{t> : (u, v ) G 7} > 0.5, e < 0.24 and 

/ s K\ /200\ 2 
(20) — >21 . 

There then exists I > [as in the definition of h in (18)/, 5\, 82, sq > suc/i 
that for G [0, do], (ao, bo) £ 7e and s G [0, so]; 

e sA ^(a /o, 60/0) > ((1 + 8 2 s)a s f s , (1 + 5 2 s)b s f s ), 

where (a s ,b s ) = (fg(s, (ao, bo)) as above, /o is defined as in (17) and / s is /o 
mi/i i/ie transition regions translated outward by 5is: 

1, x G [— L + I — 8\s, L — I + 5\s\, 

. h(x + L + 8\s), x G [-L - I - 81s, -L + I - Sis], 

{ > Js[X) ~ \ h{L + 8 lS -x), xe [L-l + 8 lS ,L + l + 8 lS }, 

0, x G (00, — L — I — 8±s) U [L + 1 + 8\s, 00). 

Remark 4.3. By abuse of notation, we again let J-„ denote the time 
evolution of the spatial system where each (u(x),v(x)), flows inde- 

pendently along the vector field n. 

For any (u,v) G TZ, the set of values {(ufo(x),vfo(x)) : x G M} forms a 
line segment in TZ with endpoints O and (u,v). Proposition 4.2 above com- 
bines the properties of the flow for small s with the spatial distribution 
generated by the heat flow. Before we prove Proposition 4.2, we state two 
technical lemmas necessary for its proof. 

Lemma 4.4. If I is fixed and f = fo is defined as in (17), then for 

x€ (-L-1- s,-L l —] U (l+—,L + 1 + s 

V 200; V 200 



and s small, we have 

e sA f(x)>f(x) + ^. 



Lemma 4.5. Let s > be fixed, fo as defined in (17) and 

i, % J fo{x) + ms, —L — l — s<x<L + l + s, 
\ 0, otherwise, 

where m > 0. Then there exist positive constants 82 depending on m but 
independent of s such that for allx, f(x) > (l + 82s)f s (x), where f s is defined 
as in (21). 
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Lemma 4.4 says that the lower part of the transition region in /o in- 
creases at a rate proportional to l// 2 , which is the case since in that part of 
the transition region, /o is convex and A/o = 0(1/Z 2 ). It can be proved 
by working with the solution to the heat equation in convolution form 
e sA g(x) = J exp(— y 2 /4s)g(x — y)dy, or working with Brownian mo- 
tions if one prefers to calculate probabilities. Lemma 4.5 involves only ele- 
mentary calculus. We refer interested readers to Yu [9] (Lemmas 5.1.6 and 
5.1.7) for full details of proofs of these two lemmas. 



Proof of Proposition 4.2. By virtue of (20), we can choose I > 
such that 

/ 200 \ 2 4.02 „ 1 



Ki > — K 2 < 



V1997 Z 2 ' 5.05Z 2 ' 
Hence, we can select m > such that 



(22) m > 



ifi/199y 2. 



2 o 01 



2 V2ooy P 

and 

(23) m -5p - 1MR 2- 

Fix 9 G [0, 8 ], (a , b ) G 70 and s G [0, Sq] for the moment. 

First, consider x G [— L — ^q,L + ^j], where the definition of /o in (17) 
implies that 

ao/o(a:) > bof (x) > \h{L + ^j) = \\{^f > 0.24 > e. 

Hence, by the first half of Assumption 4.1, 

F%(aofo(x),b fo(x)) 

> ((1 + K lS )a s f (x), (1 + K lS )bJ (x)) 

for s G [0,so] and x G [-L — ^q,L + ^q]. Furthermore, by (19), we have 
eSA fo( x ) > fo(x) — 2s// 2 for all x and s > 0. This, together with the mono- 
tonicity of e sA , shows that 

(e sA F°)(a f ,b f )(x) 

> e sA ((l + K lS )a s f (x), (1 + K lS )b s f (x)) 

> ((1 + K lS )a s (f (x) - ^) , (1 + K lS )b s (f (x) - ^)) ■ 
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Using the fact that f (x) > f(L + ^) = \{^) 2 for x G [-L - L + 
we can bound the ^-coordinate in the above as follows: 



2s 



b 8 (l + K 1S ) ( /o(x) - — )=b s [ f (x) + K 1 f (x)s - - 



2s 2K x s' 



>b f 



>b. 



' . , K X /199\ 2 2s 2K lS ? 
/o(x) + ^UoJ S ~P P- 



fo(x) + 



^1/199X2 2.01 
^~V20oj ~~P~ 



if s is sufficiently small. The ti-coordinate can be treated analogously, there- 
fore 

(24) (e sA F°)(a f ,b f )(x) > (o s (/ (z) + ms), b.(f (x) + ms)), 
by (22). 

We now consider x G [— L — I — s, —L — ^] U [L + ^q,L + 1 + s]. By the 
second half of Assumption 4.1, we have 

F°{aa , ab ) > ((1 - K 2 s)a s f (x), (1 - K 2 s)b s f (x)) 

for such x. Combining this with Lemma 4.4 and the monotonicity of e sA , 
we obtain 

(e sA ^)(a /o, hf )(x) > e sA ((l - K 2 s)a s f (x), (1 - K 2 s)b s f {x)) 



(25) 



> ((1 - K 2 s)a s (/ (s) + ^2 ) , (1 " K 2 s)b s (/ (z) + -^yj 



> a 



fo(x) + 



5/= 



l.0lK 2 )s 



fo(x) + 



5/ 2 



1MK 2 )S 



> (a s (fo(x) +ms),b s (f (x) + ms)) 



by (23), where the third inequality requires s to be sufficiently small. Com- 
bining (24), (25) and Lemma 4.5 then yields the claim. □ 

Corollary 4.6. Suppose that rj satisfies Assumption 4.1 with 
min{t> : (u, v) G 7} = D x < d\ < min{v : (u, v) G (1 + #0)7} • IfT i- s sufficiently 
large and vq{x) > D\ for x G [—L + l,L — l], then vt(x) > d\ forx G [— 3L, 3L] . 

Proof. This is an easy consequence of Proposition 4.2 and the nonlinear 
Trotter product formula (Proposition 15.5.2 from Taylor [7]) 

(ut,vt)= lim (e^ A p/ n ) n (f,g). 

n — >oo 

Here, the convergence occurs in the space i3C 1 (]R), which consists of func- 
tions / such that /' is bounded and continuous on K and both / and 
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/' can be extended to be continuous on the compactification K via the 
point at infinity; the norm used here is || • ||oo + ||9/9a;( , )|| 00 . For sufficiently 
large n, Proposition 4.2 says that application of e^- 1 ^ n ^ A J- 1 ^ n to the function 
(ao/oi ^0/0) with (a, b) G 70 yields a result that is larger than (a/, 6/), where 
(a, 6) G 7e+5/n f° r some 5 > (whose exact value depends on 62 and 0) and 
the "flat region" in / is at least 25\/n larger than that of /o- □ 

It remains to check that if c is sufficiently large, then the vector field 
r/i(u, v) = (2c(l — u) + 1)?; — u, 

(26) 

T]2(u, V) = (c(u — v) — 2)v 

as in (15) satisfies Assumption 4.1. 

Lemma 4.7. If c is sufficiently large, then we can find sq, Oq, K\, K2, 
e > and a curve 7 such that Assumption 4.1 and the assumptions of Propo- 
sition 4.2 are satisfied for the flow generated by (26 ). 

PROOF. For c > 8,800, we construct a vector field £ = (£i(t£, v),^2{u,v)) 
for (u, v) £lZ such that 

< r)i(u,v), &(u,v) <?] 2 {u,v), (u,v)eTZ 

and we will show that the flow generated by f; satisfies Assumption 4.1 
and the assumptions of Proposition 4.2 with e = 0.24, K\ = 44 and K2 = 2. 
Note that there is a fair amount of leeway in the choice of the constants, as 
we have not striven for optimality in that respect. 

Define A = (0.51,0.51), B = (0.55,0.5), C= (0.9,0.5), D = (1.0,0.5) and 
let the curve 7 be given be 7 = 70 = AB U BD (see Figure 4), where AB is 
the line segment connecting A and B. For 9 G [—0.54,0.2], let 

A e = (i + e)A, B e = (i + e)B, 

C e = (0.9, (1 + 0)0.5), D e = (1.0, (1 + 0)0.5) 

and 

70 = A e B e U BqCq. 
Hence, 70 = (1 + 0)7 n TZ. Put 

A =^4_o.54, B = 5-0.54, C = C_o.54, D' = D-0.54, 

A" = A _2, B" = Bq,2, C = C0.2, D" = Dq.2- 

If we let e' = 0.23, then i?', C" and £)' lie on the horizontal line {v = e'}, while 
B", C" and D" lie on the horizontal line {v = 0.6} and n v (A") = 0.2346 < e. 
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We define the vector field £ as a piecewise linear function in the following 
way: let F\ = 400 and F 2 = 75. This implies that 

(27) ^ T? ((i^i-8)A5.1F 2 )>if 1 , 

which will be needed later on. For (it, v) £ 1Z, let £ = (£i(tj,t;),£ 2 (7j,7;)) he 
given by 

(F lU ,-2v) mKi = {e'u<v<e',u<0.9} 
or K 2 = {e' <v< 0.6, u < l.lu}, 



[l.lF 2 u,F 2 v) on Li = 



(0, F 2 v) in ^3 = (interior of the trapezoid B'C'C'B") U C'C", 
(—u, F 2 v) in 72-4 = interior of rectangle C'D'D"C", 
(rji(u,v),r] 2 (u,v)) in the rest of 7£; 

see Figure 4. To verify that £ < rj on 72- if c > 8, 800, we first observe that 
771 > — u and 772 > — 2v are trivial bounds that hold for all of 1Z. The other 
definitions of £ can be verified region by region, as follows. 

TZi U 1Z 2 : Since v > e'u and u < 0.9, we have 

771 > (0.2c + l)e'u -u> (0.046c - 1)« > 400u, 
if c> 8, 800. 
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L\\ The validity of r/i > £1 comes from the same argument as for 
TZi U 1Z2- For £2 j since v>e' = 0.23 and u = l.lu, we have 

r]2(u, v) = (Q.lcv — 2)v > 75v, 

if c> 3, 400. 
7^.3: Since u < 0.6 and v > e'u, we have 

771 (u, t/) > (2c(0.4) + l)e'u - u, 

if c > 5. The validity of 772 > £2 comes from the same argument as 
for L\. 

IZ4: The validity of 772 > £2 comes from the same argument as for L\. 

Let s = [(log(12/ll))/F 2 ] A [(log(45/23))/(Fi + 2)]. These choices guar- 
antee that for s E [0, so], , u v {T^{uq, vq)) < 0.6 if ^o < 0.55 and 

■F|(uo,«o) = (e FlS u ,e- 2s v ) E {v > §u} 
if vo > 0.5-Uo- For 8 E [0,0.2], s E [0, so] and (uq,vo) E 70, define 

(28) (u s ,v s ) = (fe(s, (uo,vo)) = intersection of OF^Uq^Vq) and jq. 

Let AA'B'E' denote the triangle with vertices A', B' and E' = (e',e') (a 
subset of 7^2) • We claim that 

(29) a(u ,v )(£KiU AA'B'E' =► J%(auo,av )>a(l + K 1 s)(u a ,v 8 ), 

(30) a(tto, wo) e Hi U AA'B'E' Tl(au ,av ) > a(l - K 2 s)(u s ,v s ). 

Notice that this implies the required inequality in Assumption 4.1 because 
Tlx U AA'B'E' lies beneath the horizontal line v = e = 0.24. 



Proof of (29). First, consider the case (vq,vo) E AqBq. As £ is lin- 
ear in 7£i and 7^2, it suffices to restrict to 9 = and a = 1 in this case, 
that is, consider (uq,Vq) = XA + (1 — X)B E 7 for some A E [0, 1]. Note that 
(1/VT7, 4/\/l7) is a unit normal vector perpendicular to Provided 
•F|(ito, fo) is to the left of B'B", the rate of increase of Fg(uo,Vo) — (uo,vq) 
in the direction of (l/\/l7, 4/VT7) is 

(F 1 ^(^|(no,wo)),-2vr 1 ,(^|(no,wo))) ■ (1/Vl7, 4/^17)* 

(iq - 8)7r„(^|(u ,wo)) ^ (Fi - 8)0.5 ^ „ 
> = > p= > Ai, 

by (27), where we can use the lower bound 'k v {J-^{uq,Vq)) > 0.5 since the 
trajectory FHuq^vq) for (uo,Vq) E AqBq stays above AqBq. On the other 
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hand, once the trajectory J-^(uo,vq) hits B'B", the rate of increase of 
Ff(uQ,Vo) — (uq,vq) in the direction of (l/y/Vf, 4/y/Vf) is 

Tr v (Fi(uo,vo))(i.iF2,F 2 ) ■ (i/vTr,4/vTr)' 

5.1F 2 7r v (J r l(u ,v )) 5.1F 2 -0.5 Jr 
> 7= — > Ki , 



.'17 - n/17 

by (27). Therefore, 

F£(uq,Vo) > (1 + K 1 s)(u s ,v s ), 



if (it ,uo) e ^e^e- 

Now, consider (uo,vo) G BqCq, with 6* G [0,0.2] and a G (0,1] such that 
cwo > e'- Then (uq,vq) and (crao,c«;o) are in TZ^. Hence, 

jF|(ano,awo) = (cmo 5 ae F2S vo) 



until, possibly, the trajectories hit B'B" (which will happen simultaneously 
by the linearity of £ in 7^3 ), when the w-coordinate will also start to increase. 
Furthermore, for (uq,vo) G BqDq, v s = vq by the geometry of 70 and the 
definition of ip in (28), thus 

f|(au ,rao) > a(l + Kis)(u a ,v s ) 
in this case as well since F 2 > K\. Similar reasoning applies when (uq,vq) G 

!F^{auo,avo) > (auoe~ s , avoe F2S ), 

where we have equality provided ,F|(erao, avo) G IZ4. Thus, before ^(no,wo) 
hits C'C", we have (u s ,v s ) = (itoe~^ 2+1 ^ s , Vq) and 

Fj:(auo,av ) > ae F2S (u s ,v s ) > a(l + K 1 s)(u s ,v s ). 

Notice that if (auo,avo) G then the u-coordinates of ^{oluq, avo) will 
stop decreasing once C'C" is hit. □ 

PROOF of (2.5). According to the definition of £ in 72 1 U72-2, 
(31) ^{auQ,avo) = a(uoe FlS ,voe~ 2s ) 

for a(uo,fo) G 1Z\ U AA'B'E', provided the trajectories remain in 1Z\ U 
AA'B'E' , which is the case by the choice of so- Assume, first, that (uo,^o) £ 
AqBq. If ^(uo,fo) does not hit [i.e., ^(auo, avo) does not hit OB'] 

by time s, then 

^{auo,av ) = aJ T l(u ,v ), 
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hence, in particular, even 

J|(a«o,Q!Uo) > a(u s ,v s ), 



which implies (30). On the other hand, if J-^{uq,vq) does hit B'B" by time 
s, then (u s ,v s ) = (pg(s,(uo,vo)) = Bg and, hence, 

J|(au ,auo) > a(u s ,v s e~ 2s ) > a(l - K 2 s)(u s , v s ), 

which establishes (30). 

Now, consider the case (uo,vq) £ BqDq, where we always have (u s ,v s ) < 
(uq,vo) (u s decreases with s until the trajectory hits B'B" , after which both 
u s and v s remain constant). Thus, (31) implies 

J|(an ,afo) > ae~ 2s (u , v ) > a(l - K 2 s)(u s ,v s ), 

so (30) also holds in this case. The proof of Lemma 4.7 is now complete. □ 

4.2. Upper bounds: Existence of d 2 and D 2 in Condition (*). We first 
establish the following proposition. 

Proposition 4.8. // c is sufficiently large, then there exist constants 
d 2 < D 2 < 1 and T such that if Vo(x) < D 2 for x G [-L + l,L— I], then 
vt{x) < d 2 for x G [— 3L,3L] for all t>T, where (ut,vt) solves the PDE (15). 

Proof. Because of the monotonicity of (15), it suffices to consider the 
uniform initial condition uq = 1,vq = 1 and to show that Vt < d 2 for large t. 
Therefore, we need only concern ourselves with the ODE 

du 

— = (2c(l-u) + l)v-u, 

^ = (c(u-v)-2)v. 
We can bound r] 2 (u,v) = (c(u — v) — 2)v for any v > 1 — 1/c as follows: 
m (u,v) = (c(u - v) - 2)v < (c(l- U--J^ -2^jv 

= -v<-( 1- -) <0, 



if c > 1. Thus, for any d 2 satisfying 1 > d 2 > 1 — 1/c, there exists T, such 
that if ito = i>o = 1, then vt < d 2 for t>T. □ 

Proof of Theorem 1.4. If we let M = L - I, then Corollary 4.6, 
Lemma 4.7 and Proposition 4.8 show that Condition (*) from the begin- 
ning of Section 4 holds for the PDE system (15). This, in turn, implies the 
conclusion of the theorem by section 3 of Durrett and Neuhauser [3] . □ 
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